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ABSTRACT 

We present convincing evidence of unburned carbon at photospheric velocities in new observations 
of 5 Type la supernovae (SNe la) obtained by the Nearby Supernova Factory. These SNe are identified 
by examining 346 spectra from 124 SNe obtained before +2.5 d relative to maximum. Detections are 
based on the presence of relatively strong C II A6580 absorption "notches" in multiple spectra of each 
SN, aided by automated fitting with the SYNAPPS code. Four of the 5 SNe in question are otherwise 
spectroscopically unremarkable, with ions and ejection velocities typical of SNe la, but spectra of the 
fifth exhibits high- velocity (v > 20,000 km s _1 ) Si II and Ca II features. On the other hand, the 
light curve properties are preferentially grouped, strongly suggesting a connection between carbon- 
positivity and broad band light curve/color behavior: Three of the 5 have relatively narrow light 
curves but also blue colors, and a fourth may be a dust-reddened member of this family. Accounting 
for signal-to-noise and phase, we estimate that 22lg°% of SNe la exhibit spectroscopic C II signatures 
as late as —5 d with respect to maximum. We place these new objects in the context of previously 
recognized carbon-positive SNe la, and consider reasonable scenarios seeking to explain a physical 
connection between light curve properties and the presence of photospheric carbon. We also examine 
the detailed evolution of the detected carbon signatures and the surrounding wavelength regions to 
shed light on the distribution of carbon in the ejecta. Our ability to reconstruct the C II A6580 
feature in detail under the assumption of purely spherical symmetry casts doubt on a "carbon blobs" 
hypothesis, but does not rule out all asymmetric models. A low volume filling factor for carbon, 
combined with line-of-sight effects, seems unlikely to explain the scarcity of detected carbon in SNe la 
by itself. 

Subject headings: supernovae: general — supernovae: individual (SN 2005cf, SN 2005di, SN 2005cl, 
SN 2005ki, SNF20080514-002) 

1. INTRODUCTION 
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There is a well-established core population of Type la 
supernovae (SNe la ) distinguished by striking spectro- 
scopic homoge neity ( |Branch, Fisher fe Nugent|1993 Fil- 
ippenko 1997). The maximum light ultraviolet-optical- 



intrared spectra of these SNe carry signatures of predom- 
inantly intermediate-mass and iron-peak elements (often 
identified: oxygen, magnesium, silicon, sulfur, calcium, 
and iron). The consistent presence of these features in 
the spectra of the SN la core population contrasts with 
the intermittent detection of carbon lines. Such carbon 
is detected in the earliest data, but fades away by the 
time of peak brightness. 

SNe la probably arise from the thermonuclear inciner- 
atio n of carbon/oxygen white dw arfs in binary systems 
(see jHillebrandt fc Niemeyer 2000 for a review). Since 
carbon is consumed during the explosion, constraints on 
its quantity, distribution, and incidence in SN la ejecta 
help identify viable explosion mechanisms. For exam- 
ple, turbulent deflagration models predict that a large 
amount of unprocessed carbon may be left ove r (e.g., 
Gamezo et aLl|2003l [Ropke fc Hillebrandt | [2005} |Ropke| 
et ai. 2007 ) while models where the deflagration is fol- 
lowe d by a detonation phase can more completely burn i t 



(e.g., |Hoflich et al.|2002l|Kasen, Ropke fc Woosley|2009 ). 
Interpreting oxygen signatures is more ambiguous, since 
oxygen is produced as well as consumed in the explosion. 
Thus, constraints on carbon in SNe la are expected to be 
an important piece of the explosion mechanism puzzle. 

Unfortunately, multiple factors hinder an unbiased 
accounting of carbon signatures in SN la spectra. 
Challenges inherent in both discovering and coordinat- 
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ing time-series spectroscopic follow-up of early, nearby 
SNe la result in relatively poor statistics at phases op- 
timal for studying carbon (i.e., within 10 days after ex- 
plosion or before about a week prior to maximum light): 
Early observations probe the outermost layers of the SN, 
the region of the ejecta where carbon is most likely to be 
found. Optical carbon features are seldom conspicuous 
even at the earliest phases, so noise may muddle detec- 
tion in "screening" spectra of young SNe la. To achieve 
large numbers of SN candidates with a high "purity" of 
early SN la discoveries, SN surveys must cover as much 
sky as possible while maintaining a few-day cadence; 
a difficult endeavor. Such surveys with associated on- 
demand access to spectroscopy resources, like the Nea rby 



Supernova Factory (SNfactory, | Aldering et al. 2002 ) or 
the Palomar Transient Factory (PTF, |Rau et al.||2009p, 



provide a direct way to address the relative lack of spec 
troscopic time series sustained from the earliest phases. 
Such data sets provide opportunities to obtain new con- 
straints on SN la physics. By comparison, spectra ob- 
tained at and after maximum light are already plentiful 
and well-studied. 

The need to understand the SN la explosion mecha- 
nism and progenitor system is as critical as ever, espe- 
cially because of the important role SNe la play in obser- 
vational cosmology. Observations of high-redshift SNe la 
resulted in the discovery of Dark Energy, a revelation 
that has changed our view of the Universe radically over 
the past decade. The Universe is not only expanding, 
but the rate of expansion itself is accelerating, in con- 
trast to previou s expectations that it would be found to 
be decelerating (Riess et al.|1998 Perlmutter et al.|1999| 
Astier et al.||2006| |Riess et al.||2007t |Wood-Vasey et al" 



et al. 
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20071 IKowafeki et" al.||2008t |Hi"cken et al.||2009| |Sullivan 
'2011| |Suzuki et al.||2011[)l New techniques show 
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2011 Blondin et al. 2011|). How- 



ever, t. 

to stellar evolution astrophysics amplifies concerns that 
drift in emergent SN properties as a function of redshift 
could be confused with a cosmology signal. For example, 
the demographics of the underlying progenitor popula- 
tion may drift over time, leading t o a shift in measured 
SN properties such as decline rate ( Mann ucci et al. | 2005[ 



Scannapieco & Bildsten 2005 Sullivan et al.|2006[ ). Over- 
coming such problems requires building a coherent pic- 
ture of the progenitor system and explosion mechanism. 

In support of this effort, we consider the occurrence 
of unburned carbon in a large set of early phase spec- 
tra obtained by the SNfactory. O ur analysis supple- 



ments that of Parrent et al. (2011) who conducted the 
first comprehensive carbon signature survey using pub- 
licly available (and some new) SN la spectroscopy. To- 
gether, our studies find or verify several "definite" or 
"probable" C II A6580 detections in SNe la that are nei- 
ther super-Chandrasekhar candidates nor otherwise ob- 
viously peculiar. Six of 27 normal SNe la with SNfactory 
observations before — 5 d and with sufficient signal-to- 
noise (S/N) have definite C II A6580 detections. Such 
normal SNe are the foundation of the SN la Hubble dia- 
gram; studying observables correlated with the incidence 
of carbon may provide new constraints on the explosion 
mechanism, leading to new systematics controls or ruling 



out the need for some. In particular, we are interested 
in addressing the question of whether carbon-positivity 
in SNe la is idiosyncratic (distributed randomly within 
the population without regard to other observables like 
light curve shape or color) or specifically correlated with 
other parameters. For example, if the frequency (or 
strength, or persistence) of carbon signatures is somehow 
correlated with luminosity or light-curve shape, explosion 
models need to reproduce this relationship. Our results 
are the strongest indication so far of such a connection 
between carbon-positivity and broad-band photometric 
properties in normal SNe la. 

The remainder of this article is arranged as follows. 
In § [2] the initial carbon detection procedure is outlined 
and in § [3} the selected sample is presented. § H] covers 
our analysis of the pre-maximum evolution of the iden- 
tified carbon features, overall characteristics of the se- 
lected spectra, and light-curve/color properties of these 
SNe. Details on individual objects are deferred to the 
Appendix. § [5] brings together our results, provides the 
first measurement of the C II A6580 detection rate, and 
discusses our findings in broader contexts. In § [6] we 
draw conclusions. 

2. DATA SELECTION 

The spectra used in this study were obtained by the 
SNfactory between 2004 and 2009 with the SuperNov a 
Integral Field Spectrograph (SNIFS, [Lantz et al.|[2004 ). 
SNIFS is a fully integrated instrument optimized for au- 
tomated observation of point sources on a structured 
background over the full ground-based optical window 
at moderate spectral resolution. It consists of a high- 
throughput wide-band pure-le nslet integral field spectro 



graph (IFS, "a la TIGER;" |Bacon et aL||l995[ [2000 
2001[ ), a multi-filter photometric channel to image the 
field in the vicinity of the IFS for atmospheric transmis- 
sion monitoring simultaneous with spectroscopy, and an 
acquisition/guiding channel. The IFS possesses a fully- 
filled 6'.'4x 6'.' 4 spectroscopic field of view subdivided into 
a grid of 15 x 15 spatial elements, a dual-channel spec- 
trograph covering 3200-5200 A and 5100-10000 A simul- 
taneously, and an internal calibration unit (continuum 
and arc lamps). SNIFS is mounted on the south bent 
Cassegrain port of the University of Hawaii 2.2 m tele- 
scope on Mauna Kea, and is operated remotely. Spectra 
of all targets were reduced using the SNfactory's dedi- 
cated dat a reduction p i peline , similar to that presented 



in § 4 of Bacon et al. (2001). A key point is that all 
spectra analyzed and presented have been corrected for 
telluric absorption. A brief discussion of th e software 
pipeline is presented in Aldering et al. ( 2006 1 and is up- 
dated in |Scalzo et al. ( 2010[ ). A detailed description of 
the instrument and reduction pipeline will be published 
in the future. 

The starting-point for the analysis is the sample of 
SNIFS spectra obtained prior to +2.5 d with respect to 
rest-frame -B-ba nd maximum, as d etermined by SALT2 
light-curve fits (|Guy et all 120071) to SNe la observed 
on 5 or more nights. This selection yields 346 epochs 
of spectrophotometry from 124 individual SNe la. We 
have exclud ed the previously pu blished carbon-positive 
SN 2006D ( |Thomas et~aT1|2007[ ) , but will include it in 
the discussion in ^ [51 The sample includes spectroscopi - 
cally "normal" SNe la ( jBranch, Fisher fc Nugent]|1993[ ), 
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Figure 1. Identifying C II A6580 features with the help of SYNAPPS. Each panel displays a set of fits to the Si II A6355 region of 
each spectrum. Respectively, red/solid and blue/dashed curves indicate fits without and with C II opacity. The left column shows non- 
detections, the right column shows C II A6580 detection candidates. Vertical grid lines represent blueshift relative to C II A6580, and the 
shaded rectangle denotes the Na I D exclusion region. 



overlu minous SN 1991T-like and SN 1999aa-like ( |Li et al. 
2001|) S Ne la, and candida te super-Chandrasekhar-mass 



5Nela flHowell et al.|2006|). A range of light curve stretch 
values 0.65 < s < 1.16 (IGoldhaber et al.||2001|) is rep- 



resented (corresp onding post-peak decline rate 1.82 > 
Am 15 (B) > 0.73, |Guy etaL]|2007[ ). However, there are 
no subluminous SN 1991bg-hke events following the most 
strict definition (i.e., exhibiting a prominent Ti II absorp- 
tion t rough at maximum light; Branch, Fisher & Nugent 
19931. 



To help identify candidate carbon features, we system- 
atically conduct fits to the relevant portion of each spec- 



trum using the SYNAPPS code (Thomas, Nugent & Meza 
2011|). SYNAPPS uses a highly parameterized but fast 



SN spectrum synthesis technique, and as such does not 
solve the full model SN atmosph eres problem, in con- 
trast to pro grams lik e PHOENIX d Hauschil dt, Baron fc| 
|Allard||1997|), CMFGEN (|Hillier & Miller||1998p ;~or~SED0NA 
( |Kasen, Thomas fc Nugent|2(J06[ ). The other codes men- 
tioned are examples designed tor full/partial non-local 
thermodynamic equilibrium radiative transfer, enabling 
detailed post-processing of stellar explosion models for 
comparison to observations. SYNAPPS is useful in the 
opposite sense, producing constraints on explosion mod- 
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els from observations, since it identifies ion species and 
estimates their characteristic ejection velocities. Given 
the schematic radiative transfer model that SYNAPPS uses 
and its as yet unquantified systematic uncertainty, we 
employ the fit results as an adjunct in our examination 
of the spectra rather than as a direct quantitative dis- 
criminant. The radiative transfer assumptions underly- 
ing SYNAPPS match those of the well-known SYNOW SN 
spectrum synth esis program ( Fisher|2000 
& Baron 2009): Spherical symmetry, a s 



Branch, Dang 
larply-defined 

blackbody-ermtting pseudo-photosphere, pure-resonance 
line transfer under Sobolev formalism, with line opac- 
ity parameterized radially using convenient functional 
forms and in wavelength assuming Boltzmann excitation. 
SYNAPPS combines this parameterized spectrum synthesis 
calcul ation with a parallel non-linear optimization frame- 



calculatio n with a parallel non-lm car optimizat ion frame- 
work ( |Kolda|[2005| |Gray fc Kolda|[2005l IGriffin fe Kolda 



20061 to reduce the need tor tedious interactive adjust 



ment of fit parameters and to assure more systematic 
sampling of the parameter space. 

We use C II A6580 as our primary carbon indicator. We 
consider it the best option, for reasons we now describe. 
C II AA4745 and 7234 may also be useful (and h ave been 
previously identified, e.g., Thomas et al. 2007) but are 
much weaker at the temperatures ot interest, and A4745 
is very likely to be blended with other lines. Elsewhere, 
C II A4267 has been considered as a carbon tracer (e.g.. 
Branch et aT] [20031 |Howell et ah] [20061 |Thomas etahl 
20071 but it is problematic: It does not seem to pro- 
duce a reliable, distinct feature in the spectrum, and is 
vastly weaker than C II A6580 at the relevant temper- 
atures. There is some evidence for the high-excit ation 
C III A4649 line in the spectra of S N 1999aa-like (|Garj 
avini et al.|2004[ ) and SN 1991T-like ( [Parrent et al.|2(ITTf 
SNe la, but it is the only strong optical C III feature and 
is not reliably detected in the spectra of less luminous 
objects. C I lines are mostly too red or too weak (e.g., 
AA9087, 9406, 10693) for us to consider with SNIFS data 
( |Marion et al.|2006[[2009[ ). Thus, C II A6580 seems to be 
the best choice, especially considering the observational 
constraints. 



As |Parrent et al. (2011) demonstrated (their Figure 9), 
C II A6580 absorptions caused by carbon at higher ve- 
locities (much faster than that at the photosphere at the 
phases of interest) will disappear into the Si II A6355 
absorption trough where they cannot be detected reli- 
ably. Hence, inferences about carbon from the C II A6580 
notch are limited to lower velocities. Fortunately, the 
low velocity regime is what is of greatest interest for 
constraining the explosion mechanism. Carbon at the 
highest velocities may not be surprising; the presence of 
carbon is much more variable at lower velocities in exist- 
ing models. 

A final, but related, issue worthy of mention is the sys- 
tematic error that could be incurred when estimating the 
ejection velocity of carbon directly from the position of 
the absorption notch instead of wit h the aid of model - 
ing, that of center-to-limb variation ( |H5fiich et al.j l990). 
Without the aid of any kind of spectral modeling, the ve- 
locity of the carbon layer could be underestimated by up 
to a factor of 2. However, SYNAPPS (at least) handles line 
blending explicitly, and the velocity at the photosphere 
will be estimated by a fit to the full spectrum, so we can 
be reasonably confident that the carbon is located near 



or just above the photosphere. 

The wavelength range selected for initial fitting ex- 
tends from 5300 to 7300 A. The pseudo-continuum level 
and line features are simultaneously fit by SYNAPPS, and 
all fits include lines from O I, Si II, and S II. Lines from 
these ions are the most important to consider in this 
wavelength region while the SN is on the rise. Strong 
iron-peak lines only become important after the phases 
considered here. A second-order polynomial warping 
function multiplied against the synthetic spectrum is 
included to mitigate limitations of the parameterized 
model continuum, which does not affect our ability to 
detect a C II A6580 absorption "notch," a much smaller 
scale feature. Since SYNAPPS models the SN spectrum 
but not the intervening gas, we mask out rest-frame 
wavelengths between 5875 and 5915 A, a generous al- 
lowance to exclude host galaxy interstellar (or circum- 
stellar) Na I D contamination, based on the width of 
the strongest Na I D feature observ ed in the sample 
(SNF20080720-001, |SNfactory|[2008bl). All spectra pre- 



sented and fit by SYNAPPS are corrected for observer- 
frame Milky Way (Ry = 3.1) extinction using the CCM 
( Cardelli et al. 1989 ) dus t law, with color excess esti- 
mates from Schlcgel et al. ( 1998[ ). 

Three fits are performed to each spectrum to find can- 
didate C II features. The first set of fits including only 
O I, Si II, and S II provides a carbon-free baseline. The 
second set of fits adds C II with a "photospheric" Sobolev 
optical depth profile, meaning that the minimum veloc- 
ity of the profile is constrained to equal the velocity at 
the photosphere. The third set of fits substitutes a "de- 
tached" C II optical depth profile, with the minimum 
velocity of the profile unconstrained with respect to the 
photospheric velocity. In all 3 fits, the Sobolev opacity 
declines with velocity exponentially. A set of C II fits al- 
lowing Si II opacity to be detached was also performed, 
but in practically all cases SYNAPPS found a better fit 
without Si II detachment. 

Some of the targets included in the analysis have pre- 
viously been identified or suspected of bei ng carbon- 
positiv e, in particular SNF20080514-002 ( |SNfactory 
2008a). To reduce (somewhat) the risk of confirmation 



bias, the process of organizing and executing the fits in- 
cludes concealing the identity of the parent SN of each 
spectrum and its phase. This is accomplished by shuf- 
fling the spectra and assigning a unique integer identifier 
to each. The mapping from spectrum to SN and phase 
is kept in a hash file that is only consulted after the car- 
bon detections have been made. This procedure, along 
with the SYNAPPS fitting, forced us to consider more spec- 
tra more carefully. It also concealed the measured light 
curve fit parameters while identifying strong carbon fea- 
tures, ensuring that any correlation between them and 
the presence of carbon would also be uncovered without 
bias. 

For a detection to be considered definitive, we require 
that C II be detected in at least 2 spectra of an object's 
time-series, either from the same night of observation or 
across multiple consecutive nights. The final C II identi- 
fications are confirmed by visual inspection and compar- 
ison of the fits in the region of the Si II A6355 emission 
peak. Cases where the fits with and without C II are 
basically identical are considered non-detections (the ad- 
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dition of C II opacity does not improve the fit). Fits 
with C II that differ from the baseline fit, and moreover 
improve upon it, are flagged for closer scrutiny. Sev- 
eral example fits are shown in Figure m While the de- 
tailed agreement in all cases is not perfect, the SYNAPPS 
fits greatly improve our confidence in the C II identifi- 
cations we otherwise would have made purely by visual 
inspection. SYNAPPS also highlighted cases we might have 
overlooked, but in practice few of these cases ultimately 
made it into the sample, usually because they fail the 
2-spectrum cut or were less convincing. 

3. THE CARBON-POSITIVE SNe la 

The fit survey yielded 6 of 124 SNe la with definitive 
C II A6580 absorption notch detections. In one of these 
cases, SN 2005cf, there is only a single-spectrum detec- 
tion, but this SN is difficult to discard since the detection 
is corroborated by other external data sets. Of the 124 
SNe la in the sample, 4 are super-Chandrasekhar-mass 
SN la candidates (including SN 20 07if, the only clea r 
C II A6580 detection among the 4, [Scalzo et~al]|2010| ). 
Removing these leaves 5 out of 120 carbon-positive 
SNe la: SNe 2005di, 2005el, 2005ki, SNF20080514-002, 
and SN 2005cf (a special case). For our dataset, spec- 
tra having the requisite S/N generally come from SNe 
at the lowest redshifts. Due to the resulting small vol- 
ume sampled by untargeted surveys out to such redshifts, 
the sample here is dominated by SNe la from targeted 
searches. These new detections are the focus of our de- 
tailed analysis in the rest of this article. Basic target 
information is summarized in Table [l] and selected pre- 
maximum spectra of each object, up to the point where 
the carbon signatures disappear, are shown in Figure [2| 

A few characteristics of the spectra presented in Fig- 
ure [2] are worthy of specific mention. The spectra of 
SN 2005di are quite red, but there is evidence (dis- 
cussed below) that this is due to extrinsic dust extinc- 
tion. SN 2005el, SN 2005ki, and SNF20080514-002 all 
possess what appear to be pronounced flux excess blue- 
ward of about 3600 A. The spectra of these 3 objects 
are also dramatically similar in terms of detailed feature 
morphology. SN 2005cf clearly stands apart, with its 
more modest relative flux levels blueward of 3600 A, but 
broader Ca II H&K and Si II 6150 A absorptions, and 
apparent high velocity (HV) Ca II infrared triplet. 

3.1. SN 2005di 

SNIFS observed SN 2005di on 10 occasions between 
Aug. 14 and Sep. 220 The first 3 spectra of SN 2005di 
appear in Figure [2j The — 8.6 d spectrum used here is a 
co-addition of 2 consecutive exposures on the same night. 
Close examination of the time series reveals the presence 
of narrow rest-frame Na I D absorption (see Figure [3^) . 
This and the relatively red appearance of the spectral 
energy distribution lead us to conclude that SN 2005di 
is heavily extinguished by line-of-sight dust in its host 
galaxy or circumstellar environment. Since SNIFS expo- 
sure times were generally not adjusted to compensate for 
high extinction, the spectra of SN 2005di have lower S/N 
than typical SNe la observed by SNIFS at the same red- 
shift. Still, C II A6580 is detected in the first spectrum 
and possibly in the second. 

1 All dates discussed in this article are UTC. 



Initial spectroscopic classifications of S N 2005di were 
published by 2 groups. Salvo et al. (2005) reported that 
a spectrum obtained Aug. 16 resembles the spectrum of 
another SN la a f ew days past peak, and |Gal-Yam, Sand| 
& Leonard ( 2005[ ) state that their Aug. 14 spectrum sug- 
gests that SiN 2005di is "young," but neither report dis- 
cusses C II. The SNIFS quick-reduction pipeline classifi- 
cation indicated that SN 2005di was a SN la about 1 week 
before maximum, and highly extinguished with an esti- 
mated E(B — V) ~ 0.5. Revisiting our quick-reduction 
spectrum we find that the C II A6580 absorption notch 
could have been identified, but its significance was not 
obvious at the time. 

SN 2005di was used in the S Nfactory spectral f lux ra- 
tio luminosity indicator study ( Bailey et al.|2009 ). That 
study found the ratio of flux measured at 642 nm to that 
at 443 nm at maximum light (7^.542/443) to be a good pre- 
dictor of relative luminosity. SN 2005di is the object in 
that study with the largest value of 'R-s^/aaz, owing pri- 
marily to its red color. Any possible effect of C II A6580 
emission near 642 nm either as a possible source signal 
or noise in the relationship between 1Zqa2/aa3 an d Hubble 
diagram residual was not discussed, since at the time of 
publication the presence of C II had not been noticed. 
We expect our C II A6580 detection to have no impli- 
cation for the flux ratio result, since that study was re- 
stricted to spectra obtained within 2.5 d of maximum, 
when the carbon signature is gone. 

3.2. SN 2005el 

Nineteen SNIFS visits of SN 2005el were conducted be- 
tween Sep. 26 and Nov. 20. Figure [2] includes the spec- 
tra from the first 2 of these. The —6.9 d observation is 
a co-addition of the 2 observations obtained that night. 
In this case, C II A6580 is only positively detected by 
SYNAPPS on the first date of observation. However, it is 
clearly present in both high S/N spectra taken that night 
and the co-add. 
SN 2005el is a very well-observed SN (|Wood-Vasey| 



et al. 2008 Hicken et al. 2009 Contreras eTaT 2010 



Its spectroscopic normality has been remarked upon m 



)J[C 

s bl 



the literatu re previously (e.g. |Phillips et al.|2007||Hoflich 



et al. 



2010). Initial classification reports indicated that 



SN 2005el was discovered a few to several days prior to 
maximum ( |Wong et al.|2005||Morrell, Folatelli fc Hamuy| 



2005 1 but until now no remarks have been made in the 
literature regarding the presence of C II in the spec- 
tra. These classification spectra were obtained the same 
night as the first SNIFS observations. The SNIFS quick- 
reduction pipeline in use during the SN 2005el campaign 
failed to properly correct for a cosmic-ray hit coincident 
with the C II A6580 notch in one spatial element of the 2 
spectra obtained on the first night. Subsequent improve- 
ments to the pipeline later revealed the clear C II A6580 
notch in both spectra. 

3.3. SN 2005ki 

SNIFS follow-up of SN 2005ki was triggered Nov. 20. 
Six nights of observations were obtained in all, ending 
Dec. 2, around the time of maximum brightness. Follow- 
up ended prematurely due to equipment failure at the 
UH 2.2 m telescope that forced an extended offline pe- 
riod for SNIFS. Therefore, we adopt the date of max- 
imum and light-curve fit parameters from a SALT2 fit 
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Table 1 

New C II-Positive SNe la. 



Supcrnova a 


Host Galaxy b 


SALT2 Fits 


B-band Maximum 


Rest Frame Phases 




Name 






c 


Date (UTC) 


of Detection 


SN 2005di 


MGC -04-52-46 


0.025298 


-1.35 ±0.22 


0.47 ± 0.03 


2005-08-26.2 


-11.3,-8.6 


SN 2005el 


NGC 1819 


0.014910 


-2.20 ±0.18 


-0.14 ±0.03 


2005-10-03.6 


-6.9 (x2) 


SN 2005ki c 


NGC 3332 


0.019207 


-1.89 ±0.08 


-0.01 ± 0.01 


2005-12-01.2 


-10.4,-8.4 


SNF20080514-002 


UGC 8472 


0.022064 


-2.00 ±0.15 


-0.07 ± 0.02 


2008-05-26.5 


-10.0,-7.9,-5.0,-3.1? 


SN 2005cf 


MCG -01-39-003 


0.006461 


-0.38 ±0.19 


+0.07 ±0.02 


2005-06-12.9 


-9.4 



a Discovery credits (all dates UTC). Lick Observatory Supernov a Search (LOSS, [Filippenko et al.|200l(: SN 20 05di on Aug. 12.4 | |Pugh| 
|& Li|2005bl, SN 2005el on Sep. 19.5 JMadison, Baek & Li|2005l, SN 2005ki on Nov. 18.6 (Ponticello & Li|2005l, SN 2005cf on May 28.4 

' jPugn L Lifcoosa) . sn- 



SNfactory: SNF20080514-00 2~on May 14. 3 | SNfactory|2008at 
° Rcdshift credits. MGC -04-52-46 jMathewson fe Ford|l996[ . RC3.9 jde Vaucouleurs et al.|l99l) : NGC 1819 
Digital Sky Survey DR7 jAbazajian et al.|2009| ' 
c Light curve parameters from 



MCG -01-39-003. Sloan 



j NGC 3332.UGC 84 72. 
a fit to data from |Hicken et~aT\ | |2009[ and |Contreras et al-l^OlO) . 



to mo re co mplete photometry publ ished by Hickcr Tet al.| 
(2009) and |Contreras et al l pOTu] ). The first 3 spectra 
of SN 2005ki are presented in Figure [2] Relatively weak 
C II A6580 absorption is detected in the first spectrum, 
and very marginally in the second. Returning to the 
first 2 SNIFS quick-reduction spectra of this object, the 
presence of C II A6580 now seems obvious in retrospect. 
Our classification announcement made no mention of it 



(SNfactory 1120051) 



3.4. SNF20080514-002 

Twenty nights of spectroscopy of SNF20080514-002 
were obtained with SNIFS, starting May 16 and ending 
Jul. 9. The first 4 spectra of SNF20080514-002 appear in 
Figure|2] The —7.9 d and —5.0 d spectra are co-additions 
of 3 consecutive observations on each respective night. 
C II A6580 is clearly detected in the first 2 spectra of 
this object, and probably in the third. Though the emis- 
sion peak visible at the position of the notch in the fourth 
spectrum looks flattened, we hesitate to posit a flux de- 
pression due to C II opacity in this case. At this point, 
SYNAPPS is able to convincingly fit the feature with or 
without a weak C II A6580 line. 



As initially reported ( SNfactory 2008a ) , the early spec- 
tra of SNF200805 14-002 are quite blue, indicative of a 
SN la at or before a week prior to maximum light. Both 
spectra exhibit absorption lines attributed to C II A6580, 
but C II A7234 is not as readily apparent as in, say, 
SN 2005el. C II A4745 was also mentioned, but is a 
much weaker and more blended line than the others. 

Like SN 2005di, SNF20080514-002 was used in the 
Bailey et al. ( 2009[ ) flux ratio study. In marked con- 
trast to SJN 2005di, this target has an extremely blue 
color but is well within the core of the Hubble residual 
versus 7^42/443 distribution (their Figure 2). As with 
SN 2005di, the use of the at-maximum spectrum for the 
flux ratio study limits the influence of C II A6580 on 
those results. 

3.5. SN 2005c! 

SNIFS observed SN 2005cf on 12 nights from Jun. 3 to 
Jul. 8. The first 2 spectra are shown at the bottom of 
Figure [2] A small C II A6580 absorption notch is visible 
in the —9.4 d spectrum but becomes absent or nearly 
absent 2 days later. Our selection guidelines require de- 
tection in 2 distinct SNIFS spectra, so we should tech- 
nically disqualify SN 2005cf from further consideration. 
However, SN 2005cf is an object with extensive previ- 



ously publish ed spectroscopy ( Garavini et al.j2007 |Wang 
et al. 2009b[). Neither previous spectrosc opic study nor 



the o riginal classification announcement ( |Modjaz et al 
2005 1 discussed the putative C II A6580 notch visi ble in 
the spectra. Wh en considering those observations, Par- 
rent et al. (2011) placed SN 2005cf into their "uncertain" 
C II detection category. But given 3 independent data 
sets showing the same feature at roughly the same phase, 
it now seems difficult to overlook the presence of C II in 
the case of SN 2005cf. 

4. ANALYSIS 
4.1. Spectroscopy 

The detailed evolution of the carbon features in the 
pre-maximum spectra of SNe 2005di, 2005el, 2005ki, and 
SNF20080514-002 is depicted in FiguregJ SN 2005cf and 
a non-detection case for comparison (SNF20071021-000) 
appear in Figure |4j The C II A6580 notch is readily ap- 
parent, situated at or quite near the rest-frame emission 
peak of the Si II A6355 feature. The detailed evolution 
of the feature in each case is described in the Appendix. 

In the first 4 cases, the notch is strongest in the earliest 
spectrum. It weakens gradually with time, and without 
appreciable wavelength shift. Generally, the notch is cen- 
tered on an apparent blueshift of 12,000 km s _1 relative 
to rest-frame C II A6580 (also roughly the same veloc- 
ity as observed in SN 2006D). In SN 2005cf, the notch 
disappears between the first and second spectrum. 

Excluding SN 2005cf, we see that while the notch dis- 
sipates, the shape of the spectrum just redward evolves 
as well. The slope of the decline in flux on the red side 
of the Si II A6355 emission feature becomes steeper while 
the notch weakens and the underlying Si II A6355 emis- 
sion peak rises. In some cases (SNe 2005di, 2005ki, and 
SNF20080514-002) the earlier spectra may even exhibit a 
mild "flux plateau" morphology. Along with an absorp- 
tion notch, such an emission plateau could be a signa- 
ture of a "detached" distribution of C II line opac ity in a 
shell separated s omewhat from the photosphere (Jeffery 



& Branch 1990). However, under spherical symmetry, 
such a plateau should be roughly symmetric about the 
line rest wavelength. Visually, this does not seem to be 
the case here. 

The notch in the spectrum of SN 2005cf is morpho- 
logically distinct from the other cases, and appears to 
span a slightly smaller range in blueshift (10,000 to 
13,000 km s _1 ). When the notch disappears, there is at 
most only a small amount of evolution in the flux slope 



Carbon in SNe la 7 



i — ■ — ■ — ■ — i i i i i 




3000 4000 5000 6000 7000 8000 1000°0 
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Figure 2. Pre- maximum SNIF S spectra of carbon-positive SNe la. Spectra up to the phase when C II A6580 dissipates are included. A 
smoothing filter ( Savitzky & Golay 1964 ) applied to the flux values is overlaid on the spectra of SN 2005di (in red) for presentation. Phase 
labels denote rest-trame phase with respect to f?-band maximum. Spectra are labeled by SN name to the right. 
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■(a) SN2005di 
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(d) SNF200805 14-002 
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Figure 3. Spectroscopic evolution of the 5800-7500 A region of the carbon-positive SNe la 2005di (a). 2005el (b), 2005ki (c), and 
SNF20080514-002 (d). Smoothing-filtered flux values are overlaid (red line) for presentation purposes. Blue vertical grid lines denote 
blueshift relative to 6580 and 7234 A, the rest wavelengths of the 2 strongest optical C II lines. The Na I D region is highlighted for 
SN 2005di, where narrow rest-frame line-of-sight absorption is clearly visible. 
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Figure 4. Probable C II detection case SN 2005cf (a) and ex- 
ample non-detection case SNF20071021-000 (b). Blue vertical grid 
lines denote blueshifts relative to C II lines as in Figure [3] In (b) 
the red overlaid curve is the smoothed flux. 



just red ward. 

Confident identification of C II A7234 is problem- 
atic. In almost none of the cases presented is there 
a clear absorption feature at the same blueshift as the 
C II A6580 notch. At most, there is only an inflection vis- 
ible around 6900-7000 A and a slight emission "bump" 
around 7100 A. Visual inspection alone would seem to 
be insufficient to draw more firm conclusions about the 
presence or absence of C II A7234 in general. 

To examine the extent to which the selected targets are 
spectroscopically normal SNe la, we again use SYNAPPS, 
this time to fit the full wavelength range of each pre- 
maximum spectrum of each selected SN (detailed com- 
mentary is provided in the Appendix). Close examina- 
tion of the C II features in the full fits may also enable 
some explosion-model independent statements about the 
distribution of carbon in the SN ejecta. Two selected fits 
for each SN appear in Figures [5] and [6] One fit to the first 
available spectrum, and another fit to the first spectrum 
where C II is absent. In general the opacity profiles for 
C II lines decrease with time as expected. Wavelengths 
of major ion contributors to the spectra are marked at 
the top of either figure. Except for the presence of C II, 
the spectra of the first 4 carbon-positive SNe la are any- 
thing but extraordinary — typical SN la ions at normal 
pre-maximum ejection velocities. SN 2005cf is another 
matter, since HV Si II and Ca II compone nts are manda- 
tory for a good fit, as disc ussed previously ( Garavini et al. 
2007| |Wang et al.||2009b[ ). 

In each case, we find compelling evidence for opacity 
from C II, O I, Mg II, Si II and Si III, S II, Ca II, singly- 
ionized ions of Ti through Co, and Fe III. For both Si II 
and Ca II we enable 2-component fits: One for a photo- 
spheric (undetached) component, and another for a HV 
(detached) component. The opacity profiles of all other 
ions are undetached from the photosphere. Optimal val- 
ues for the velocity at the photosphere are found in the 
range 12,000 to 13,000 km s -1 . Including some Na I 
opacity does seem to improve the fit morphology around 
the Si II 5800 A absorption feature — however this iden - 
tification is not unprecedented (e.g., Branch et al.|2 005). 
In all cases, we were able to reconstruct the basic struc- 
ture redward of 3500 A quite well, and at times had rea- 
sonable success fitting the ultraviolet (UV, where lines 
from iron-peak elements dominate). We do not expect 
exceptional agreement in the UV, as the underlying as- 
sumptions of SYNAPPS (in particular the treatment of 
line opacity) begin to break down here. In general, the 
strengths of the features evolve smoothly with time, and 
no abrupt changes are observed as the C II becomes op- 
tically thin. 

SYNAPPS accounts for the C II A6580 feature, includ- 
ing the region just redward of the notch, rather well (see 
Figure Yn red curve) . Variations with and without de- 
tached TJ II were tried, but again SYNAPPS tended to set 
the lower velocity of the C II region to the velocity at 
the photosphere, making the need for detachment un- 
compelling. Analysis of the fit residuals indicates that 
the SYNAPPS model has insufficient fidelity for achieving 
a reduced \ 2 of 1, but we find no strong indications of 
persistent residual structure. 

In Figure [7J the 5800-7500 A region of the first spec- 
trum of each carbon-positive SN is overlaid with the con- 
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Figure 5. SYNAPPS fits (red lines) to selected pre-maximum spectra of SNe 2005di, 2005el, and 2005ki (panels a, b, and c respectively). 
The fit to the first available spectrum of each SN is shown, followed by the first fit where C II opacity is no longer invoked. Major 
contributions to the fit are labeled by ion in panel (a). Singly ionized atoms of Ti through Co contribute to the spectrum blueward of 
3500 A. 
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Figure 6. SYNAPPS fits (red lines) to selected pre-maximum spectra of SNF20080514-002 (a) and SN 2005cf (b). The fit to the first 
available spectrum of each SN is shown, followed by the first fit where C II opacity is no longer invoked. Major contributions to the fit are 
labeled by ion in panel (a). Singly ionized atoms of Ti through Co contribute to the spectrum blueward of 3500 A. 
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Figure 7. Action of C II opacity in converged SYNAPPS fits to 
the earliest observed spectrum of each SN. The red/solid curve 
plots the synthesized spectrum from the best fit with all ions. The 
blue/dashed curve results when the final synthetic spectrum is re- 
computed without C II. The carbon opacity screens the emission 
peak from Si II, in addition to creating the notch shape. 



verged, all-ions, full-wavelength SYNAPPS fits (shown in 
red). The blue curve results when the fully converged 
synthetic spectrum is recomputed with the C II opac- 
ity deactivated, while all other settings retain their con- 
verged values. The blue curves are thus not re- fits with- 
out C II included, in contrast to the example plots shown 
in Figure IT] Thus, the change from the blue curve to the 
red one illustrates what role C II plays in the formation 
of the spectrum quite clearly. Not only is C II A6580 re- 
sponsible for the absorption notch, it "screens" emission 
from bluer Si II A6355 effectively enough that some emis- 
sion relative to the C II-free spectrum (centered around 
6600 A) is visible in the top 4 cases. Finally, Figure W 
indicates that the flux depression around 6900-7000 A 
may indeed be influenced by C II A7234 to some extent. 

4.2. Light Curve and Color Properties 

It is interesting that 3 of the 5 carbon-positive SNe la 
have relatively narrow light curves, with low SALT2 x\ 
values around —2 (s ~ 0.87, or Ami5(B) ~ 1.47, see 
Table fl]) . These same targets also have bluer than nor- 
mal colors. SNF200805 14- 002 is also among the SNe la 



observed by Swift/UVOT (|Gehrels e t al. 2004; R oming 
et al. 2005[) with largest UV excesses ( jlmmler &~brown 
2008 1 , as shown in Figure [8] (adapted from Milne et al. 
50TU{. Swift/UVOT photometry for SNe 2008hv and 



5NF200805 14-002 will be presented in Brown et al. 2011, 
in preparation. Figure|9]is a plot of SALT2 fit parameters 
for all SNe la with SNIFS light curve fits and a spectrum 
before —5 d, the latest phase for which any C II A6580 
is positively detected in this sample. SNe 2005el, 2005ki, 
and SNF20080514-002 form a relatively tight "cluster" at 
the fast-declining (low xi) but blue corner of this plot. 
Other published res ults for SN 2005el have a somewhat 
larger value for x\ (|Hicken et al 



2009 Contreras et al. 



20101, though all measurements indicate that X\ is neg- 



ative and this SN belongs in this group of relatively nar 
row light curves with blue colors. Thus, the quantitative 
results for SN 2005el may be further refined, but our 
qualitative conclusions should remain unchanged. 

Neither SN 2005cf nor SN 2005di fit into the fast- 
decliner/blue-color family. In the case of SN 2005cf, this 
may not be surprising since the Si II and Ca II features 
quite robustly trace out HV ejecta more dramatically 
than in the other SNe at the same phases, indicating 
real physical differences between SN 2005cf and the oth- 
ers. However, as described in § [3] there is evidence that 
SN 2005di is extrinsically reddened. Since our data do 
not resolve the Voigt profiles of the possible multiple ve- 
locity components comprising the Na I D absorption line, 
we are unable to apply a quantitative c orrection for ex- 



tincti on using Na I D equivalent width (Poznanski et al. 
20111. Instead, we use a different procedure to inves 



tigate whether SN 2005di could simply be a reddened 
member of the fast-decliner/blue-color subclass. 

To see if this is at least a plausible hypothesis, we 
warp the spectral time series of SN 2005di according to 
the CCM dust law, fitting for the optimal differential 
E(B — V) that brings the spectra into detailed agree- 
ment above 5000 A with each of SNe 2005el, 2005ki and 
SNF200805 14-002. We make no attempt to constrain the 
differential Ry and we assume a standard Milky Way 
value of 3.1. If we assume that the main difference be- 
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Figure 8. SN la colors from Swift/UVOT. Note that SNF20080514-002 (open hexagon), SN 2008Q (right filled semicircle), and SN 2008hv 
(left filled semicircle) are among the objects with the largest UV excess observed by Swift. In particular, the objects are very clearly separated 
in w — v color space from the rest of the sample. 

3 cases is shown. The agreement above 5000 A on a 
feature-by-feature basis is quite good, though there is a 
significant difference in the Si II A6355 absorption fea- 
ture with respect to SN 2005ki — clearly an "intrinsic" 
difference. Blueward of the 5000 A cut, the extrapo- 
lated color looks quite consistent. Thus, while we can- 
not prove definitively that SN 2005di is a dust-estranged 
member of the fast-decliner/blue-color family, there does 
seem to be a resemblance. Figure [9] includes the position 
of SN 2005di when corrected for the median derived dif- 
ferential E(B — V) = 0.40 is removed. 

5. DISCUSSION 
To summarize, our most important findings are: 

• Definitive new detections of C II A6580 absorp- 
tion notches in the pre-maximum spectra of 5 
largely spectroscopically normal SNe la, at appar- 
ent blueshifts around 12,000 km s _1 , consistent 
with the velocity of the ejecta at the photosphere. 

• The observed notches fade with time without ap- 
preciable shift in wavelength, disappearing be- 
fore maximum light. Detailed fitting shows that 
C II A6580 influences the formation of the spec- 



tween SN 2005di and the other SNe is the amount of 
dust extinction the former suffers, then correcting for 
it should bring the spectral time series into alignment. 
In the case of perfect "twin" SNe, the agreement above 
5000 A should be exquisite on a feature-by- feature basis, 
and the extrapolated flux below 5000 A a good match. 

The warping technique and a fully dev eloped twin SN 
study are the subject of a future paper ( jFakhouri et al. 
2011 1, here we briefly describe the mechanics tor 2 given 
objects. Spectral observations of the objects are paired 
based on phase, with similar phases being matched to 
each other. A linear interpolation correction to the fluxes 
accounts for overall flux difference arising from the small 
phase differences. A CCM dust law with fixed Rv = 3.1 
is fit to the ensemble of paired spectra, as well as an indi- 
vidual scale factor for each pair, to account for any small 
scatter. By fitting all of the pairs (extending from early 
epochs to late epochs) jointly for the appropriate dust 
law warping, we minimize the influence of any particular 
features on t he fit at any given phase. 

Figure [10] shows pairs of spectra that result when the 
time series of SN 2005di is warped to match the flux 
from another SN above 5000 A, roughly matching by 
phase. The earliest pair of phases formed in each of the 
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Figure 9. SALT2 fit parameters for SNe la followed with SNIFS having observations earlier than —5 d with respect to maximum light 
(narrower SN la light curves correspond to negative xi, broader ones to positive x\). Red open squares mark the positions of carbon-positive 
SNe la, blue circles mark the others. The green open triangle represents SN 2005di after removing a differential E(B — V) obtained by 
warping the spectral energy distribution to match those of SN 2005el, 2005ki, and SNF20080514-002. Three of the 5 identified carbon- 
positive SNe la form a family of fast-decliner/blue-color objects. Other SNe la (labeled by name, but blue circles) with similar show very 
marginal evidence for C II A6580, but their first phase of observation is around —7 d. 

trum beyond the immediate vicinity of the absorp- 
tion notch. 

• Three of the 6 carbon-positive SNe la (including 
SN 2006D) in the full SNfactory sample group to- 
gether in light curve shape and color space at the 
fast-decliner/blue-color corner. A further SN may 
belong to this family, plausibly separated from it 
by the effect of extrinsic dust reddening. 

We consider 4 main questions in light of our results. 
First, factoring in the limitations of our sample and ap- 
proach, what is the incidence of photospheric carbon in 
pre-maximum spectra of normal SNe la? Second, to 
what extent does carbon-positivity in SNe la predict a 
fast-declining light curves and blue color (or vice- versa), 
and what is the physical significance of such a connec- 
tion? Third, what inferences can we make about the 
geometry of the ejected carbon in the SNe la studied 
here? Finally, what impact do our findings have on our 
understanding of the SN la explosion mechanism? 

5.1. Incidence of C II 

Our data suggest that unburned carbon, while far from 
ubiquitous, is also not particularly rare — in line with 



the result of Parrent et al. (2011 ). A detailed rate calcu- 
lation is beyond the scope of this article, but we can try 
to estimate the frequency of photospheric-velocity car- 
bon in normal SNe la based on the latest phase (—5.0 d) 
and lowest median S/N per 2.4 A bin in the wavelength 
range of interest (11.7, that of SN 2005di at -11.3 d) 
where C II A6580 is detected. 31 SNe la pass these cuts 
in the full SNfactor y sample (which includes SN 2006D, 
Thomas et al. 2007), of which 2 are super-Chandrasekhar 
candidates and a further 2 are SN 1991T-like. Excluding 
these 4 SNe places the fraction of SNe la with C II A6580 
that persists to —5 d at 22+g % (Bayesian binomial con- 
Ca~meron||2010[ ). 



fidence interval, see 



5.2. Carbon in SNe la 

SN 2005el, SN 2005ki, SNF20080514-002, and perhaps 
SN 2005di could be members of a physically meaningful 
family of SNe la characterized by detectable carbon in 
their early spectra, relatively fast-declining light curves, 
and intrinsically blue colors. These 4 SNe were selected 
solely based on the presence of carbon signatures in their 
spectra, so we now reverse the selection and re-examine 
the early spectra of other objects with similar light curve 
or color parameters. 

The most interesting results from our surveyed sam- 
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Figure 10. Selected spectra extracted from time series of SN 2005di, warped by a CCM dust law to match the time series of SNe 2005cl, 
2005ki, and SNF20080514-002. SN 2005di at -8.6 d is compared to SN 2005el at -6.9 d at top; at -11.3 d is compared to SN 2005ki at 
-10.4 d in the middle; and at -11.3 d to SNF20080514-002 at -10.0 d. The red line is the wavelength portion of the SN 2005di time-series 
(above 5000 A) where the data were warped to match. The blue line is SN 2005di, extrapolated with the derived CCM dust law correction. 
The black line is a spectrum from the target SN warped to in each case. 

inary decline-rate estimate from the UVOT photometry 
appears to be faster than average. 

Next, we consider whether normal SNe la with previ- 
ous C II A6580 detections in the literature have a ten- 
dency to wards narrow l ight c urves and blue colors. A 
revisit of Parrent et al. (2011 1 yields mixed results, par- 
tially because of the lack of consistent light-curve fits 
across all the data sets. Of the 2 "definite" cases pre- 
sented there, SN 1990N fails to fit the pattern (xi = 
0.54 ± 0.10, c = 0.10 ±0.01); the other is SN 2006D. The 
"probable" category (where data are available) is domi- 
nated by narrower light curves than average as well, but 
no strong tendency toward blue colors. The e xact role 
that host ga laxy dust extinction plays in the |Parrent| 



shows the 5800-7500 A region of the pre- 
maximum spectra of SN 2006dm (Schwehr & Li 20061 
and SNF20080919-002. Neither case presents a robust 
notch, however the Si II A6355 emission peak looks some- 
what "squashed" in the earliest spectra. The evidence is 
hardly decisively in favor of photospheric-velocity car- 
bon in these 2 additional S Ne la, but it is tant alizing. 
Carbon-positive SN 2006D ([Thomas et al.||2007[ ) is also 
plotted in Figure [9j It also nas a SAL1'2 xl around —2 
but is redder than the objects uncovered by our survey 
here. 

Returning to the Swift /UVOT c olor-cu rves in Figure 8l 
we identify SN 2008Q ( [Villi et al.|[2008| and SN 2008hv 
(Pignata et al. 2008) as having U V colors similar to 
those of SNF20(Ml5T4-002. All 3 of these targets sep- 
arate from the rest of the Swift sample especially well 
in u — v color space. While the spectroscop ic classifi- 
cation announcements f or neither SN 2008Q (Stauishev 
nor SN 2008hv flMarion, Garnavich fe GerarcTy 



et al 
all 



2008 



2008 
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mention C II A6580, the pre-maximum spectra of both 
objects include the absorption notch signature (V. Stan- 
ishev and G. H. Marion, priy ate c ommunication) . Both 



Ganeshalingam et al. (20101 and Brown et al. (20101 
report a taster-than-average decline rate (correspond- 
ing roughly to x\ of —0.92 and —1.67 respectively) for 
SN 2008Q, suggesting that it could also belong with the 
fast-decliner/blue-color family. Optical light curve anal- 
ysis for SN 2008hv is yet to be published, but a prelim- 



(2011) sample is unconstrained, however. Over- 
would seem that the picture for unburned carbon 
in SNe la may be a complex one, where there are multiple 
subclasses and considerable variation between them. 

Four of the 5 newly identified carbon-positive SNe la 
may indeed be representative of a meaningful subclass, 
but SN 2005cf stands apart. Its light curve shape is less 
extreme and its colors more red. It also lacks the large 
UV excess of SNF20080514-002 seen in Figure [8} HV 
ejecta signatures in the spectra of SN 2005cf are much 
stronger than in the other carbon-positive SNe la con- 
sidered here. Though HV Ca II features may be "ubiq- 



uitous" (Mazzali et al. 2005), it is clear that there is 
considerable variation in the velocities and line strengths 
observed. Cases with the strongest HV Ca II IR triplet 
also seem to exhibit strong Si II A6355 absorption early 
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Figure 11. C II A6580 feature regions for SN 2006dm and SNF20080919-002, 2 SNe la with pre-maximum coverage and light curve 
widths similar to those of SNe 2005el, 2005ki, and SNF20080514-002; however these objects are redder. Distinct C II A6580 notches are 
not obviously visible in these objects, but the apparent Si II A6355 emission peaks appear flattened early on. 



Oil UStanishev et al.|2007| . The origin of the HV features 
remains unclear (Is it a density or composition enhance- 
ment in the ej ecta, or is it swept-up circumstellar mate- 
rial? See e.g. JGerardy et al.| [2004l |Tanaka et aT1[2008l 
Wang et al. 2009b). Evidence for a connection between 
carbon and the HV features is weak. The best candidate 
aside from SN 2005cf with HV features and a possible 



C II A6580 notch seems t o be SN 2003du (Stanishev et al. 
2007[|Tanaka et al.|2008| ). On the other hand, SN 2009ig' 



(also with HV features), a ppears to be a counter -example 



(not a distinct notch, so |Parrent et al. 
at most a 



2011 label it 
SNF2007102T-000, the 



"possible" detection 
non-detection example in Figure [4] is yet another 

A number of "peculiar" SNe la (or SN la-like objects) 
with carbon lines have been discussed in the literature. 
Here we compare our sample of normal SNe la with those 
objects. 

C II features are a spectroscopic bellwether of candi- 
date super-Chandrasckhar SNe la, where they seem to 



persist to maximum light and even beyond (|Howell et al 
20061 jHicken et al.||2007| IScalzo et ai.||201(J[ [ Yamanaka 
et al'p009[ [Silverman et al.||2010| ). This has been inter- 
preted as evidence supporting the super-Chandrasekhar 
hypothesis, where more total carbon is present in the 
progenitor and/or much remains after disruption. These 
objects are quite overluminous, and have much broader 
than normal light curves. Their spectra, while mostly 
consistent with one another, are quite different from 
those of objects considered here. Overall, if the super- 



Chandrasekhar candidates form a subclass of SNe la, it 
is distinct from the one proposed here. 

With varying degrees of certainty, carbon is also oc- 
casionally detected in low-lumino sity, short time-scale, 
or oth erwise peculiar SNe la. |Taubenberger et al.| 
(20081 claim that C II A6580 is present in the spectra 
oTSN 2005bl, a low-luminosity SN 1991bg-like SN la, 
but carbon is seldom remarked upon or observed in 
other members of this subclass. More unusual events 



20031 iBranch et al. 
2006;|Philiips et al. 



like SN 200 2cx jLi et al.. 
SN 2005hk dChornock et al' 
SN 2007qd (|McClellancreraT.P7TD|)', and SN 2008" 



ley et al. 2009' 



|2010aj [Valenti et al. 



2001 , 

Sum 



2009) have been 



conjectured to represent a distinct population of pecu- 
liar, low-luminosity SNe la. These SNe are quite sub- 
luminous (peak Mb > — 18), have light curves that 
evolve more rapidly than normal SNe la, and have spec- 
tra that exhibit many narrow absorption lines that sug- 
gest lower kinetic energies than typical. Pure white dwarf 
deflagrations or failed deflagrations may provide a plau- 
sible way to account for those observations, in which 



case carbon signatures may also be expected (e.g., Ropke 
2007| ). The carbon identifications in the spectra of 



et al. 



these events are irregular, tentative in most cases, though 
SN 2008ha pres ents the most clear C II signatures (Fo 



ley et al.||2010a|) . Spectra of candidate "SNe la" (Bl 
st'en et al. p'UUTp SN 2002bj (|Po znanski et al.||2010r"ah7i 



SN 24J10X ( |Kasliwal et aT| 2010| ) exhibit robust C hde_ 
tections, but so tar detailed simulations of SNe .la ( Shen 
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et al-pOlO ) do not predict observable carbon. All of the 
above cases are much dimmer than the SNe considered 
here, and muc h more exotic spect roscopically. 



SN 2006bt ([Foley et al 
cific 



2010b), however, merits spe- 
An outlier on the SN la light-curve- 
SN 2006bt exhibits some 



attention, 
width/luminosity relation, 
spectroscopic peculiarities. Most interestingly, its ear- 
liest spectra possess a unique C II A6580 signature, a 
notch centered at 5,200 km s _1 , apparent in at least 
2 consecutive observations. There is little, if any, ev- 
idence for emission to the red of the notch, and the 
vast difference between the apparent blueshift of the 
C II A6580 notch and the velocities of other ejected ele- 
ments (12,500 km s " 1 for Si II A6355) is quite striking. 
Foley et al. (2010b I suggest that the discrepancy can be 
accounted for by an ejected blob moving at an angle to 
the line of sight. SN 2006bt at the present seems to be 
the strongest evidence for carbon ejecta distributed non- 
spherically in SNe la. 

5.3. Geometry 

Inspired by the characteristics of multi-dimensional ex- 
plosion simulations with strong asymmetries, observers 
have searched for spectroscopic and spectropolarimct- 
ric indications of deviation from spherical composition 
symmetry in SNe la. Spectropolarimetry offers the best 
means for probing the geometry of SN ejecta, but oppor- 
tunities for follo w-up are relatively rare an d the available 
sample is small. Wang & Wheeler ( 2008 ) describe spec- 
tropolarimetric observations of 2 SINe la in our carbon- 
positive sample, SNe 2005cf and 2005el (a single epoch 
each). In neither case is there a particularly strong case 
for deviation from basically spherical composition sym- 
metry. A dominant axis in the Stokes Q/U diagram is at 
least identifiable for SN 2005el, suggestive of some com- 
bination of global and smaller scale asymmetries. 

While spectropolarimetry may be more definitive for 
constraining the geometry of the unburned material, it is 
interesting to consider the constraints we can set with our 
pure flux measurements. The comparison of synthetic 
spectra with and without C II shown in Figure [7] suggests 
that the spatial distribution of unburned carbon in the 
ejecta may be quite consistent with spherical symmetry. 
Figure 12 demonstrates our reasoning simply, using toy 
model synthetic line profiles. 

The classic SN spectroscopic feature formed by a single 
line is plotted at the top of Figure[l2] Blucshifted absorp- 
tion produced by material in front of the photosphere and 
expanding toward the observer, with an emission hump 
centered at the rest wavelength of the line arising from 
material adjacent to the photosphere in projection. For 
discussion purposes, we label this line "line A" and assign 
it a rest wavelength of 6355 A. 

Consider a second line, "line B" with rest wavelength 
6580 A, a Doppler shift of 12,000 km s _1 , a shift quite 
similar to the cases seen in the data. If the ejecta bear- 
ing the parent ions of line B are distributed in a "blob" 
that only partially covers the photosphere and/or line A 
emissivity in projection, it is quite easy to construct a 
line B absorption notch blended with the line A emission 
component. An examp le o f the resulting blend appears 
in the middle of Figure 12 overlaid on the line A profile. 

Now, suppose the line B ions are distributed spherically 




Line A only 
(spherical symmetry) 

Line A with 
nonspherical Line B 

Line A with 
spherically symmetric 
Line B 



6000 6400 

Wavelength (A) 



6800 



Figure 12. Simulated, blended line profiles for different distribu- 
tions of opacity and line strength. The black dotted curve traces a 
baseline spectrum, arising from a single unblended line ("line A") 
at rest wavelength 6355 A. In all cases, line A opacity is spher- 
ically symmetric. The red/solid curve is the spectrum resulting 
when line A is blended with another line ("line B" 6580 A) and 
line B opacity is distributed in a blob only partially blocking the 
photosphere. The blue/dashed curve shows the blend with line B 
opacity distributed in spherical symmetry, and demonstrates the 
screening effect observed in the carbon-positive SNe la discussed 
here. The red/solid and blue/dashed curves are overlaid on the 
unblended line A profile for visual reference, as was done with the 
synthetic spectra with and without C II in Figure [7] 

symmetrically. In this quite different blended 

profile results. Line B notches can be produced, but not 
without additional disruption of the line A profile. Flux 
redward of the notch is suppressed and there is even a 
mild red emission feature relative to the single unblended 
line A case. Such a blend is depicted at the bottom of 



Figure 12 again overlaid on the single line A profile for 
comparison. 

Because line B opacity does not encircle the photo- 
sphere in the blob case, its effect is restricted to wave- 
lengths corresponding roughly to where the blob is opti- 
cally thick in blueshift parallel to the line of sight. How- 
ever, when line B ions are distributed under spherical 
symmetry, the photosphere is cocooned by them, and 
they scatter away continuum and line A photons alike 
over a larger range of projected velocity. This results 
in the notch, the flux suppression to the red of it, and 
the small emission further on. Essentially, the driving 
factor is the projected area over which line B photons 
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screen continuum but also line A emissivity, at a given 
wavelength. 

The fits in Figure [7] appear more consistent with the 
spherically symmetric case demonstrated in Figure |12| 
Our toy models are quite simple, but illustrate the 
leading-order screening effect quite clearly. At the very 
least, our results highlight the critical need for detailed 
spectroscopy and spectropolarimetry of SNe la at early 
phases to make the most definitive conclusions about the 
geometry of the carbon ejecta. 

5.4. Implications for Explosion Models 

The possible existence of a family of SNe la with 
similar properties to those demonstrated by SN 2005el, 
SN 2005ki, SNF20080514-002, and perhaps SN 2005di, 
is particularly intriguing. The combination of a narrow 
light curve, large UV excess or at least blue colors, and 
the presence of photospheric velocity carbon suggests a 
very particular conceptual picture. A faster than average 
decline rate indicates lower than average luminosity for 
a SN la, in turn suggesting a lower than average mass 
of freshly synthesized radioactive nickel. The blue colors 
suggest that the amount of iron-peak elements synthe- 
sized or mixed into the outer layers is low (suppressing 
line-blanketing at shorter wavelengths), and these ele- 
ments are limited to lower velocities. Also, a relatively 
high carbon abundance in the material above the pho- 
tosphere indicates less efficient processing of the white 
dwarf progenitor material. 

The above picture is somewhat consistent with a subset 
or region of delayed-detonation model space. Estimates 
of the mass fraction of photospheri c carbon in SNe la 



are generally low (of order 1%, e.g. Branch et al. 2003 
Marion et ah] [20061 | Thomas et aLl|2007[ |Tanaka et al". 



2008 Maeda et al. 



2010a), and are also generally con- 



sistent with most such models. Of particular interest 
is the recent grid of multi-dimension al models presented 
by Kasen, Ropke & Woosley (2009). In these models, a 
strong ignition (arising from numerous initiating sparks) 
raises the amount of burning that takes place during 
the initial sub-sonic deflagration phase. This allows the 
outer, unburned layers to pre-expand more relative to the 
weaker ignition case. When the later detonation wave 
passes through the pre-expanded material, nuclear burn- 
ing is less efficient and the abundance of iron-peak ele- 
ments is lower relative to the weak ignition case. This 
results in relatively lower radioactive nickel mass which 
is more centrally peaked in the ejecta. 

However, the total amount of unburned carbon left in 
these models is below 1% by mass, and the S/N of the 
synthesized Monte Carlo spectra is too low to attempt 
to detect any carbon signatures in them. This result is 
actually at variance to e arlier (one-dimensio nal) delayed- 
detonation models (e.g., |Hoflich et al.|2002[ ), where some 
amount of unburned carbon can be left behind. Simply 
put, the multi-dimensional model grids may not yet be 
extensive enough. 

In most of the carbon-positive SNe la discussed here, 
there seems to be intriguing evidence that the carbon is 
distributed roughly spherically symmetrically, and at ve- 
locities overlapping silicon. There may still be room for 
deviation from spherical symmetry in the distribution of 
unburned material that we cannot rule out here. For 
example, large scale asymmetries, such as proposed by 



Maeda et al. (2010b) to explain spectroscopic diversity, 
may allow tor almost hemispherical covering of the pho- 
tosphere by carbon. In that particular picture, an offset 
ignition in the white dwarf progenitor yields an asym- 
metric distribution of burning products. This scenario is 
proposed to explain the origin of "high- velocity gradient" 
(HVG) and "low-velocity gradient" (LVG) subclasses of 
SNe la, segregated by the rate of decrease in measured 
Si II A6355 blue-shift with time ( |Benetti et aLl|2005| ). 
Viewed from the side closest to the ignition, one ob- 
serves an LVG SN la; and viewed from the other an HVG 
SN la. Detonation products are distributed in the outer 
layers in all directions, but it is no t yet clear how com- 
plete the burning is. Interestingly, Tanaka et al. (2008) 
remarked that of 3 LVG SNe la they studied, all had 
marginally detected C II A6580, and blue colors. Their 
analysis suggests that LVG SNe la undergo less intense 
burning than their counterparts. A systematic exami- 
nation of Si II A6355 blueshift evolution of SNfactory 
SNe la is in progress, but our preliminary results indicate 
that all 5 objects identified here are al so LVG members. 
Whether or not a Maeda et al. (2010b) scenario is capa- 
ble of explaining all the observables, and can yield carbon 
amounts and distributions consistent with observations, 
remains to be seen through detailed multi-dimensional 
explosion modeling and radiative transfer. 

6. CONCLUSION 

We have presented new detections of carbon signatures 
in the spectra of several SNe la. For the most part, the 
carbon-positive SNe are normal spectroscopically, with 
the only significant deviation being the presence of HV 
Ca II and Si II in one of them. Four of the 5 objects dis- 
cussed may form a robust subfamily of SNe la with a rel- 
atively fast decline rate and intrinsically blue colors. The 
carbon is detected at velocities consistent with the ejec- 
tion velocity at the photosphere, and in the same velocity 
range as other freshly synthesized material. The distri- 
bution of carbon is apparently consistent with spherical 
or near-spherical symmetry. The behavior of most of 
the carbon-positive SNe la discussed here suggests strong 
radial abundance stratification, consistent with some de- 
layed detonation models. 

Only one of the objects uncovered in the course of our 
analysis has ever been previously identified in the litera- 
ture as carbon-positive, either at the time of initial spec- 
troscopic classification or in subsequent analysis. Our 
intent in pointing this out is not to criticize; indeed we 
have highlighted cases that earlier escaped our notice as 
well. Rather, we suggest that the C II A6580 notch sig- 
nature is not well-known, easy to overlook, or easy to 
mistake for a reduction artifact. In fact, this is part 
of why we required 2 spectra for labeling detections as 
definitive. The examples provided here should help to 
sensitize observers to carbon signatures in early-phase 
SN la spectra. 

Understanding the frequency, distribution, and 
amount of unburned carbon in SNe la is indeed a use- 
ful discriminator among competing explosion models, so 
constraints on theory depend on reliable and efficient 
early triggering of spectroscopic (and spectropolarimet- 
ric) follow-up of SNe la for detailed study. As current 
SN la cosmology efforts approach their limits in terms 
of calibration and empirical standardization, enhancing 
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our understanding of SN la physics becomes all the more 
important. Wide-field and high-throughput robotic sur- 
veys have the power to grow the early-time data set, 
when coupled efficiently with ample spectroscopic follow- 
up resources. This opens new opportunities for fully ex- 
ploiting the analysis of the outer layers of SN la ejecta 
to constrain explosion models. 
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APPENDIX 



NOTES ON INDIVIDUAL SUPERNOVAE 
SN 2005di 

C II. The C II A6580 absorption notch is clearly visible in the —11.3 d spectrum, located between roughly 9,000 and 
15,000 km s" 1 with a minimum around 12,000 km s _1 in apparent blueshift. Almost 3 days later, the notch is barely 
discernible, and it is clearly absent by —3.6 d. In the last spectrum, it is unlikely that the noise is obscuring a notch: 
The change in the shape of the spectrum between the last 2 spectra is overt. Examining the C II A7234 region, we 
detect an inflection in the shape of all 3 spectra, but its position is inconsistent across phase. This inconsistency may 
be a consequence of the low S/N of the observations. In the first spectrum, the flux depression is coincident with the 
A6580 notch in blueshift, but this is not the case in the later spectra. However, in all 3 spectra, a bump in the flux 
centered at a A7234 blueshift of 5,000 km s _1 is clearly visible. Figure [3^, also clearly shows the narrow rest-frame 
Na I D absorption, in the shaded rectangle. 

SYNAPPS Fits. C II contributes to the spectrum on —11.3 d and also in the next spectrum on —8.6 d, but on 
-3.6 d SYNAPPS finds no need for it. The SYNAPPS velocity at the photosphere (v phot ) fits are 12,600, 12,200, and 
11,900 km s _1 respectively. HV Ca II opacity is invoked (at v > 20,000 km s _1 ) though the signature is rather small 
it is more important for fitting the IR triplet than the UV H&K feature. Ultimately, the case for HV Si II is 
uncompelling, which is not surprising since the blue side of the Si II A6355 absorption only extends to 18,0 00 km s , 
in contrast to instances where HV Si II has been previously discussed (e.g., SN 2005cf, Garavini et al.|20"07 ) and found 
to extend upwards of 20,000 km s _1 . 



SN 2005el 

C II. A C II A6580 absorption notch is readily seen in the —6.9 d spectrum. As in SN 2005di, the notch is located in 
the apparent blueshift interval 10,000 to 15,000 km s _1 . There is a clear depression in the flux around 6980 A, possibly 
attributable to C II A7234 absorption, though the blueshift is not completely consistent with that for the A6580 notch. 
Both features disappear from the spectrum by the time of maximum light. As in SN 2005di, the red side of the Si II 
emission feature steepens as the SN ages and the C II A6580 dissipates. There is also an emission bump to the red 
of the A7234 notch, centered at a rest-frame blueshift of 5,000 km s _1 . Like the 2 notches, this bump disappears by 
maximum light. 

SYNAPPS Fits. SYNAPPS requires C II only in the —6.9 d spectrum, and not at maximum light. The value of v p hot 
changes from 12,500 to 12,000 km s _1 from the first spectrum to the second. HV Ca II and Si II opacity components 
were enabled as was done with SN 2005di. The results were similar, with some HV Ca II contributing to the IR triplet 
and the bluer feature of the H&K absorption blend, and HV Si II making only a minor contribution to the blue wing 
of the Si II A6355 absorption in the first spectrum. 
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SN 2005ki 

C II. In the —10.4 d spectrum, the C II A6580 notch appears slightly more blueshifted than in the first spectra 
of SNe 2005di and 2005el, extending just beyond 15,000 km s _1 . The same is true for Si II A6355 absorptions. As 
before, the carbon signatures dissipate as the SN brightens: Weaker by —8.5 d and absent by —5.5 d. Also, in the first 
spectrum, an emission plateau just redward of the notch is visible, asymmetric with respect to the rest wavelength of 
A6580, and again it weakens in parallel with the notch. A clear association for C II A7234 in this spectrum is quite 
problematic: The fading bump at 7100 A is again visible, but there is no clear absorption notch. 

SYNAPPS Fits. C II opacity is used in the -10.4 and -8.5 d spectra, but not in the -5.5 d spectrum. SYNAPPS 
photospheric velocities for SN 2005ki are a bit higher than in the other SNe considered here, with v p hot fit at 13,100, 
12,600, and 12,500 km s _1 . One major difference is that a distinct HV Ca II opacity component is not required to 
reproduce the Ca II features, especially at —10.4 d: A single opacity profile extending to high velocity is sufficient. 

SNF20080514-002 

C II. The C II A6580 notch (apparent blueshift interval again 10,000 to 15,000 km s _1 ) is most apparent at —10.0 
and —7.9 d, still discernable at —5.0 d, and probably absent at —3.1 d. The notch is accompanied by a modest and 
fading emission plateau to the red which is almost symmetric about the line rest wavelength. The A7234 blueshift 
region is again characterized by an inflection and bump centered blueward of 7234 A, and again the strength of this 
bump decreases with time. 

SYNAPPS Fits. C II opacity is needed in the -10.0, -7.9, and -5.0 d fits, but not on day -3.1. The values of 
Vphot fit are 12, 700, 12, 400, 12, 100 and 12, 100 km s _1 , respectively. Only a modest amount of HV Ca II opacity is 
used in these fits, and HV Si II proves to be unnecessary. 

SN 2005cf 

C II. The C II A6580 absorption notch shape is significantly different from that seen in the other carbon-positive 
SNe la studied here. The notch is more narrow, extending in apparent blueshift from 10,000 to 13,000 km s _1 in the 
—9.4 d spectrum. The notch as faded away by —7.4 d, a phase at which it is still present in the other carbon-positive 
SNe studied here. 

SYNAPPS Fits. SYNAPPS uses C II for the —9.4 d spectrum but not for the spectrum obtained a week later. 
SYNAPPS determines v p hot = 12,800 and then 12,000 km s _1 , respectively. The —9.4 d spectrum requires both HV 
Si II and Ca II, though the HV Si II feature becomes optically thin by —2.4 d. The HV component is also clearly 
quite important for the formation of the Ca II H&K absorption in both spectra. The strength of the HV features in 
SN 2005cf distinguishes it from the other SNe la considered here. 
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